Anisotropy, Itineracy, and Magnetic Frustration in High-Tc Iron Pnictides 
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Using first-principle density functional theory calculations combined with insight from a tight- 
binding representation, dynamical mean field theory, and linear response theory, we have exten- 
sively investigated the electronic structures and magnetic interactions of nine ferropnictides rep- 
resenting three different structural classes. The calculated magnetic interactions are found to be 
short-range, and the nearest (Jia) and next-nearest (J2) exchange constants follow the universal 
trend of Jia/2J2 ~ 1, despite their itinerant origin and extreme sensitivity to the z-position of As. 
These results bear on the discussion of itineracy versus magnetic frustration as the key factor in 
stabilizing the superconducting ground state. The calculated spin wave dispersions show strong 
magnetic anisotropy in the Fe plane, in contrast to cuprates. 
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Recent discovery of the new high-temperature super- 
conductor, LaOi_a;Fa;FeAs with a transition temperature 
(Tc) of 26X [H has triggered tremendous research activ- 
ities on iron pnictides. Rare-earth {RE) doping increases 
Tc up to b?>K for Sm [2,, ji]. Replacing RE-0 layers with 
Li produces an intrinsic superconductor LiFeAs with Tc 
of 18K The 122 ferropnictides, ^LFe2As2 (AL: Ca, 
Sr, Ba, KV, span another structural class with Tc up to 
38K [I,0,l3,[l,|^,[l3]- More recently, arsenic- free FeSei_ 5 
and Fe(Sei_2:,Tea;)i_5 without any interlayer between Fe- 
(Se,Te) planes were found to be superconducting at Tc 
as high as 27-?^ under pressure [Hi [H, Q • In spite of 
the accumulating reports of both experiments and theo- 
ries, the nature of the superconductivity and magnetism 
is still far from clear. After several works have ruled out 
the electron-phonon coupling [l^, and the coexis- 
tence of magnetic fluctuation and superconductivity be- 
ing confirmed by //SR [17], intensive investigations have 
been focused on the magnetic properties of these sys- 



tributions, as well as their anisotropy and the proximity 
of the ratio Jia/2J2 to unity, which has been discussed 



terns 18 



25| . From the studies up 
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to now, one of the common and evident features is the 
interplay between superconductivity and magnetism. It 
is clear, from the different structures, that the essential 
physics lies in the iron plane forming the 2-dimensional 
spin lattice. 



Except for the Fe(Se,Te) family suggested to have dif- 
ferent magnetic structures by recent studies 18, 19, 20], 
it is widely believed that the first three classes of Fe pnic- 
tides have a common superconducting mechanism closely 
related to magnetic interactions. In order to clarify the 
raised issues and lead to further understanding, it is of 
key importance to investigate the exchange interactions 
across different classes of compounds and examine any 
trend or common features. However, material-specific 
information of magnetic interactions is scarce in spite 
of active research efforts. The direct probe of spin dy- 
namics is inelastic neutron scattering, which has been 
recently performed for SrFe2As2 [9] and CaFe2As2 [1^. 
They have revealed that the combination of nearest and 
next nearest neighbor exchange interactions \Jia + 2J2] 
is about lOOmeV, but detailed data from individual con- 



21, 22, \M, are stiU 



extensively in recent publications 
missing. 

In this Letter, using first-principle linear response cal- 
culations 



27 



28| , we provide the data of in-plane mag- 
netic exchange couplings for several Fe-based supercon- 
ductors, and discuss their spin wave dispersions. The 
data bear on the question of whether the values of ex- 
change constants indicates magnetic fluctuations play an 
important role. A total of nine materials have been 
studied: RE FeAsO {RE: La, Cc, Pr, Nd), ^LFe2As2 
{AL : Ca, Sr, Ba, K), and LiFeAs. Exchange interac- 
tions of these systems are found to be short-range de- 
spite the metallic density-of-states (DOS), and the calcu- 
lated interaction strengths follow the universal behavior 
of Jia ~ 2 J2 for all materials, a relation that arises inde- 



pendently in the frustrated magnetic picture [21|, |22| . 
Considering not only the variety of the materials studied 
here but also the high sensitivity of the Fe moment to the 
z-position of As atom 2^ 3^ , this universal behavior of 
the exchange interactions is impressive. The calculated 
spin-wave dispersion shows an anisotropic spin interac- 
tion which is different from the cuprates. 

There have been several published tight-binding (TB) 
parametrizations of the electronic structure of prototyp- 
ical LaOFeAs in the vicinity of the Fermi level using fits 
based either on Wannier functions or atomic basis sets 
33L 3^ . However the current situation still looks 



31, 32 



complicated because the projected DOS deduced from 
electronic structure calculations are based on the spher- 
ical harmonic projectors within the atomic spheres that 
may not be very well suited for the extended Fe and As 
orbitals presented here. Due to these complications even 
the crystal field splitting of Fe d level appears to be con- 
troversial in the current literature 31 , 33, H, 34 ) . 

To better understand the complicated electronic struc- 
ture around Fermi level, we performed TB analysis by 
considering dxz and dyz orbitals of Fe t2g manifold hy- 
bridizing with the arsenic Px and Py, respectively. As 
shown in Fig. [U the separation between the energy levels 
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FIG. 1: (Color online) The tight-binding band structure of 
LaFeAsO. The circles at the Fermi level on the F and M 
points indicate the hole and electron pockets, respectively. 



of Fe-t2g and As-px.y states is about 1.6eV. Account- 
ing for the hybridization matrix element between d^z- 
Px , dyz-Py states, which is of the order of l.SeV, pro- 
duces bonding and antibonding bands, both having the 
bandwidth of 2.8eV with the Fermi level falling into the 
antibonding part of the spectrum (approximately leV 
above the Fe t2g level). We also take into account the 
dxy state of Fe which hybridizes with itself (hopping in- 
tegral is approximately O.SeV), which produces an ad- 
ditional bandwidth of 2.2ey. The resulting bandwidth 
of Fe d-electron character near the Fermi level becomes 
2.8 -I- 2.2/2 = 3.9eV as exactly seen in the LDA cal- 
35| . The coordinate system used for this TB 



culation 

description is the original crystallographic lattice where 
the spin alternates in the (7r,7r) direction. In this pic- 
ture, the F-centered hole pockets (small circle in Fig. [T]) 
are mostly of dxy character, and the M-centered pockets 
(large circle) are of dxz, dyz character. This picture can 
be fine-tuned further by including the d^^-yi state which 
lies O.SeV below the Fermi level and hybridizes primar- 
ily with As-px,y states (hopping integral is about 0.8eV) 
as well as hybridization between dxz,yz orbitals with As 
Pz states (hopping integral is about QAeV). Note that 
in this picture the Fe d^a.]^ orbital becomes unoccupied 
and lies \eV above the Fermi level. 

Now we discuss the exchange interactions. To cal- 
culate the interactions between Fe moments, we used 
linear response theory [s^ [s^] based on first-principle 
density functional theory (DFT) calculations, which has 
been successfully applied to the 3d transition-metal ox- 
ides and the 5/ actinides metallic alloys [s^,!!^. We used 
the full potential linearized muffin-tin orbital (LMTO) as 
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FIG. 2: (Color online) Spin arrangment and exchange inter- 
actions in the Fe plane of the striped Qm-AFM phase. The 
arrows on lattice sites indicate the Fe spin directions. 



the basis set [39| and local spin density approximation 
(LSDA) for the exchange-correlation (XC) energy func- 
tional. The LSDA is fairly good to describe the itinerant 
Fe Sd states in these materials as shown in the previous 
studies and comparisons with angle resolved photoemis- 
sion [H H m, lioj. In the calculations of REOYeAs 
compounds, we used the LSDA-hDMFT method [23, HI] 



in which the RE 4/ orbitals are treated as the localized 
ones within Hubbard I approximation. U = QeV and 
Jh = 0.86ey were used as the on-site Coulomb repulsion 
and Hund's rule exchange parameter. Lattice constants 
are taken from experiments, and we performed the calcu- 
lations at various z(As), including experimental z(yls)oxp 
and LDA optimized z{As)lda- 

Fig. [2] shows the spin structure of the Fe plane which 
is common to the all these materials. From here on we 
use the (tt, 0) striped AFM coordinate system, which is 
convenient to discuss the spin wave dispersions. Mag- 
netic interactions between Fe moments are governed by 
two dominating AFM couplings Jiq and J2, and the FM 
nearest-neighbor exchange Ju is small. We found the 
exchange couplings J(q) can be expressed in terms of 
short-range exchange constants. This character suggests 
pursuing a comparison with local moment models with 
AFM spin interactions 21, 2^, 2^. The short range cou- 



plings do not conflict with the itinerant magnet picture 
because although the Fe 3d orbital has flnite DOS at the 
Fermi level, the magnetic interactions can still remain 
short range, which possibly reflects the bad metallicity 
and some correlation effects. 

The calculated Fe magnetic moments and exchange in- 
teractions are summarized in Table HI We use the con- 
vention that positive J means AFM couplings. The cal- 
culated moments are consistent throughout the materi- 
als. The calculations done at experimental z(As)exp are 
known to predict the moments about twice as large as 
experimental values, while at optimized z{As)lda they 
give smaller moments. The cases in which DFT overesti- 
mates magnetic moments are rare, and the cause is still 
under debate for Fe oxypnictides. Although some theo- 
rists suggest it is due to the frustrated magnetic struc- 
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System 


Moment 


Jla 


J2 


Jib 


Jla/2J2 Jla + 2J2 


LaFeAsO 


1.69 


47.4 


22.4 


-6.9 


1.06 


92.2 


CeFeAsO 


1.79 


31.6 


15.4 


2.0 


1.03 


62.4 


PrFeAsO 


1.76 


57.2 


18.2 


3.4 


1.57 


93.6 


NdFeAsO 


1.49 


42.1 


15.2 


-1.7 


1.38 


72.5 


CaFe2As2 


1.51 


36.6 


19.4 


-2.8 


0.95 


75.4 


SrFeaAsa 


1.69 


42.0 


16.0 


2.6 


1.31 


74.0 


BaFe2As2 


1.68 


43.0 


14.3 


-3.1 


1.51 


71.5 


KFe2As2 


1.58 


42.5 


15.0 


-2.9 


1.42 


72.5 


LiFeAs 


1.69 


43.4 


22.9 


-2.5 


0.95 


89.2 



LaOFeAs 

1.4 -SrFeAs 



TABLE I: Calculated Fe moments (in fig) and in-plane ex- 
change interactions (in meV), using experimental z{As). 



ture 21 1 , Mazin and Johannes suggest an alternative pic- 
ture 4l| based on magnetic fluctuation and inhomogeni- 
eties. Importantly, the electronic structure features such 
as electron-hole symmetry and the exchange interaction 
strengths are better described with z{As)exv when com- 
pared to available experimental data [2^ . Thus our 
discussion will be based on the results from z{As)exp- 
The sensitivity of moments and exchange interactions to 
z{As) is large. For example, in LaFeAsO the change of 
z{As) by DMA {Az{As) = 0.005 in terms of internal co- 
ordinates) induces about 10% difference in the moment 
and up to 20% in the exchange interactions [2^. The 
same order of sensitivity was also reported for CaFe2As2 
[sot - Therefore the deviation of up to 8% for moments 
and 30% for major exchange interactions ( Jia and J2) are 
not significant, and become much smaller if z{As) could 
be refined for each material. Taking this into account, we 
can say that the magnetic moments and exchange inter- 
actions are uniform throughout the materials considered 
here. 

One of the most important quantities to understand 
the magnetism and possibly the superconducting mech- 
anism in these materials is the ratio of Jia/2J2 , which 
has so far not been measured nor calculated. According 
to the spin Hamiltonian models [2l|, assuming Fe 

pnictides as magnetic Mott insulators like cuprates, at 
Jia/2J2 ~ 1 the system is close to the quantum critical 
regime, so a superconducting ground state may appea r 
as a result of the magnetic fluctuation [2l|, [23. |42| . 
Note that the calculated ratios shown in Table |T] are all 
around unity, demonstrating that this universal behavior 
of Jia/2 J2 can arise from itinerant magnetism, without 
the system being close to a Mott transition. The devi- 
ations of t/ia/2J2 from unity reflect not only the intrin- 
sic material properties but also the sensitive dependence 
on z{As). Although there is no apparent relation be- 
tween the Jia/2 J2 ratio and Tc, the universal feature of 
Jia/2J2 near unity is closely associated to superconduc- 
tivity since it is present throughout the materials studied 
here. The connection between itinerant AFM and super- 
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FIG. 3: value (exp.) and theoretically optimized value (opt.) 
of z{As). In the inset figure, the LDA-optimized z(As) is set 
to be zero reference. 



conductivity has been discussed previously [43|. 

Another important quantity is |Jia + 2J2| which de- 
termines the spin wave velocity in the (7r,0) direction, 
and can be directly probed by neutron scattering exper- 
iments. The available experimental data are in general 
agreement with our calculation. For SrFe2As2 calcula- 
tion shows I Jla -f 2J2I = 74meV, not much smaller than 
the 100 ± 20meV measured by neutron scattering 
Also, for CaFe2As2 our calculated \Jia + 2J2I = 75meV 
is slightly smaller than the measured 95±16meV [2^ (de- 
rived from the observed spin wave velocity, see eq. ([2]) be- 
low). Especially our result for BaFe2As2 is in good agree- 
ment with recent experiment by Ewings et al. [44] : One 
of their best fits shows that Jia — 36meV, J2 = 18mel^, 
and Jif, = —7meV. 

As an example. Fig. [3] shows the z(As)-dependence of 
the magnetic moments and interactions of SrFe2As2. The 
moment is a simple monatomic function of z{As) rang- 
ing from 0.35/ig to 2.23fj,g. The three J's have different 
behaviors. Jia increases rapidly with z{As) at the be- 
ginning, saturates in the middle, and eventually turns 
down. This behavior is a result of the hybridization be- 
tween Fe dxz,yz and As Px,y orbitals, as we discussed in 
the tight-biding representation. Due to the shape and 
orientation of the Fe dxz,yz and As Px.y orbitals, there 
is a certain z{As) that gives the maximum overlapping, 
and hence largest Jia- Note that Ju changes sign at 
z(^s)oxp, and eventually surpasses J2- Also, Jia and 
J2 plateau in the small region around z{As)exp- Sim- 
ilar behaviors are also found in other materials. The 
Jia/2J2 and \Jia + 2J2I values presented in Table U are 
robust against the small deviations in z{As) around the 
experimental values. From the data one can also calcu- 
late Jia/2J2 vs. z{As), which reveals the existence of 
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FIG. 4: The calculated spin wave dispersion of SrFe2 As2 along 
high-symmetry lines, the exchagne constants are given in Ta- 
ble H 

the "sweet spot" where the optimal ratio Jia/2J2 = 1 is 
achieved independent of any Heisenberg model assump- 
tion. In the case of SrFe2As2 it is z{As) = 0.357 (the 
inset of Fig.E]). 

The calculated spin wave dispersion gives more in- 
tuitive information about the magnetic interaction and 
anisotropy of these systems [2^ . The dispersion rela- 
tion of the 2D striped- AFM lattice reads 

c^(q) = S^iJo + Jib(q))2 - (./i„(q) + J2(q))2. (1) 

Using the calculated magnetic exchange constants, we 
plot the spin wave dispersion of SrFe2As2 in Fig. lU 
whose S — 0.94 is taken from experiment [8]. The non- 
symmetric dispersions in (0, 0) — (0, tt) and (0, 0) — (tt, 0) 
directions indicate in-plane magnetic anisotropy, which is 
a major difference from cuprates. At small q near (0,0), 
the spin wave velocity in the (7r,0) direction is 

v± = 2aS \Jia + 2J2\ , (2) 

which is the relation used to experimentally determine 
[Jia + 2J2|, such as for SrFe2As2 0. The difference 
in Jia and Jib, a direct consequence of the Qm-AFM 
ordering that breaks in-plane symmetry, accounts for 
the anisotropy in (7r,0) and (0,7r) directions. These 
anisotropic spin waves can be directly probed by neutron 
scattering experiments on single crystals. 

To conclude, we have studied magnetic exchange in- 
teractions in various Fe-based high Tq superconductors 
using first-principle based linear response calculations. 
From the nine different materials, the magnetic inter- 
actions are short-range and can be well described by the 
first and second nearest-neighboring interactions. Impor- 
tantly J1/2J2 is close to unity for all the cases, just as 
would be the case for the frustration limit of a local mo- 
ment model. Calculated spin wave dispersions show the 
magnetic anisotropy and the roles of the three in-plane 
exchange interactions. Our result strongly suggests the 
magnetic fluctuation as the pairing mechanism for the 
superconducting ground state. 
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